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Abstract 

« ; 

Scalar meson contributions in chiral quark condensate are calculated in the ana- 
lytically regularized Nambu - Jona-Lasinio model using the framework of mean-field 
expansion in bilocal-source formalism. The sigma-meson contribution for physical val- 
ues of the parameters is found to be small. Pion contribution is found to be significant 
! and should be taken into account for the choice of the parameter values. 

m 1 Introduction 



Nambu-Jona-Lasinio (NJL) model [1] was the first field theoretical model of dynam- 
ical chiral symmetry breaking (DCSB) in hadron physics. The NJL model has been 
reformulated in quark language [2] during the seventies and eighties, and since then 
Qh! it exists as a successful effective model of quantum chromodynamics of light hadrons 

in the non-perturbative region. Subsequently, the NJL model was intensively inves- 
tigated also as a model of hadron matter at finite temperature and density [3]. (For 
more references see also reviews [4] and [5].) 

In overwhelming majority of these investigations, the NJL model has been consid- 
^ ■ ered in the mean- field approximation (Hartree approximation), or in leading order of 

l/n c -expansion ( n c is the number of colors). The successes in the phenomenological 
applications have led to an analysis of the structure of the NJL model in next-to- 
leading order of the expansions ( see [6]-[10] and refs. therein). Such an analysis is 
necessary for clarification of the region of applicability of results and the stability under 
a variation of parameters and quantum fluctuations due to higher-order effects. 

Since the NJL model in the mean-field approximation includes quark loops, the 
essential aspect of application of this model is a regularization. Most common regular- 
izations for NJL model traditionally entail a four-dimensional cutoff in Euclidean mo- 
menta or a three-dimensional momentum cutoff. Other regularization schemes (Pauli- 
Villars regularization or non-local Gauss formfactors) also are used for the NJL model. 
The least common regularization for NJL model is a dimensional regularization (Thus 
in reviews cited above, [4] and [5] this regularization is not even mentioned). At first 
it might appear somewhat strange, since the merits of dimensional regularization are 
generally known, and the same regularization most commonly used for calculations in 
renormalized theories, particularly in gauge theories. Apparently this fact is connected 



* Physical Department of Baku State University, Baku, Azerbaijan 

^Institute for High Energy Physics, Protvino, Moscow region, Russia (e-mail: rochev@mx.ihep.su) 







with the following circumstance: in contrast to renormalized models, a parameter of 
regularization in the NJL model is included in formulae for physical quantities, and 
it is one of the essential parameters of the model. But the parameter of dimensional 
regularization, which is traditionally treated as a deviation in physical dimension of 
space, does not permit any physical interpretation in this treatment. 

However, an alternative treatment of dimensional regularization exists - as a vari- 
ant of an analytical regularization. In this treatment all calculations are made in four- 
dimensional Euclidean momentum space, and the regularization parameter is treated 
as a power of a weight function, which regularizes divergent integrals. Such treatment 
of dimensional regularization, based on ideas of Wilson and Collins, is consistently de- 
veloped and applied to the NJL model by Krewald and Nakayama [11] in the mean- field 
approximation. It should be stressed that in this treatment of dimensional regulariza- 
tion, the regularization parameter is not at all a deviation in the physical dimension 
of space. 

We suppose that a possible treatment of this parameter is a power of some factor, 
which is a measure of gluon influence on the effective four-fermion quark self-action 
of NJL model. In some sense this possible treatment has something in common with 
non-local variants of the NJL model (see [10], [12]). 

In this work we study the NJL model with dimensional regularization in the treat- 
ment of Krewald and Nakayama 1 in the next-to-leading order of the mean-field ex- 
pansion. To formulate the mean-field expansion (Section 1) we have used an iteration 
scheme of solution of Schwinger-Dyson equation with fermion bilocal source, which has 
been developed in the literature [13]. Analytical regularization of the NJL model is 
discussed in Section 2. A purpose of our calculation is to study meson contributions to 
quark chiral condensate, which is a principal order parameter in models of DCSB. As 
our calculations demonstrate (Section 3), the pion contribution to chiral condensate 
is expressed in the analytical regularization by a very simple formula (29)- a ratio of 
pion contribution to leading approximation condensate is inversely proportional to the 
regularization parameter and does not depend on other parameters of the model. The 
pion contribution is rather significant, and tends to infinity at a bound of admissible 
values of the parameter , i.e., the model is unstable with respect to quantum fluctu- 
ations near this bound. The sigma- meson contribution is small for admissible values 
of the parameter. In Section 4 the results of Sections 2 and 3, which were obtained 
for the classical variant of the NJL model with U(l)-symmetry, are generalized for a 
physically interesting model with two flavors and n c colors (SU(2)-model). 

A choice of physical values for SU(2)-model parameters is discussed in Section 5. 
A discussion of the results is presented in the Conclusion. Quite briefly our result can 
be stated as follows: for physical values of parameters, the analytically regularized 
NJL model does not contain any pathological quantum fluctuations, connected with 
the scalar-meson contributions, though the pion contribution is significant and should 
be taken into account in phenomenological treatments of NJL-type models. 

lr To avoid unnecessary associations with the usual treatment of dimensional regularization we shall in 
what follows, refer to this regularization as analytical regularization of NJL model 
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2 Mean-Field Expansion in Bilocal-Source Ap- 
proach 

By the U(l)-model, we mean a theory of self interacting spinor field ip defined by the 
Lagrangian 

£ = ^ + |^V) 2 + (^75V') 2 )- (1) 

Here g > is a coupling. This Lagrangian is invariant under transformations of the 
chiral group Uy{l) x Ua(1)- Such a model has no direct physical applications, but, as 
is shown in Section 4, the results of physical SU(2)-model for meson contributions in 
chiral condensate are practically identical to the results of the U(l)-model. This is not 
surprising since these contributions are purely dynamical. Distinctions are displayed 
as simple coefficients (see Eq. (40) of Section 4) . 

A generating functional of Green functions (vacuum expectation values of T-products 
of fields) can be represented as the functional integral with bilocal source 

G{rj) = f D(ip,ip)ex-pi{ J dxC - J dxdytp(y)r](y,x)t(;(x)^ , (2) 
where rj(y, x) is the bilocal source of the spinor field. 

The functional derivative of G over source rj is a one-particle (two-point) Green 
function (a propagator of the field tp)'- 

() ° = i<0\TU(x)^(y)}\0>=S(x-y). (3) 

??=o L J 



5r](y,x) 

The n-th functional derivative of G over source r/ is the n-particle (2n-point) Green 
function: 



5 n G 



Sv(yi,xi)---5rj(y n ,x n ) 



T]=0 



| t{v(xiMj/i) • • • i>{x n )4>(y n )} | >= S n 



xi yi 



n 



Translational invariance of the functional-integration measure in Eq. (2) give us the 
functional-differential Schwinger-Dyson equation (SDE) for the generating functional 
G: 



5{x-y)G + id x — - + igl— rtr— 7-75^-7 V tr75 r7 \ \ = ( A > 

oi]{y,x) ^dr]{y,x) dr]{x,x) or](y,x) oij(x,x)j 



= / dxirj(x, xi)- 
J ( 



5G 



Srf(y,xi) 

We shall solve this equation employing the method proposed in [13]. (See also the brief 
review in [14].) 

For the NJL model this method is a variant of the mean-field expansion. 
A leading approximation is an approximation of the functional-differential SDE (4) 
by the equation with zero r.h.s: 

^ x 5r](y,x) ^^5r](y,x) Srj(x,x) ^ 5 5r](y,x) ^ 5 5r](x,x)i 

A solution of Eq.(5) is the functional 

G(°) =exp{Tr(^°)*r / )} (6) 
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(here Tr is a trace in operator sense, and * is operator multiplication), where is a 
solution of equation 

5(x) +idS(°\x) +ig{s (0 \x)tTS^(0) - j 5 S^\x)tTj 5 S (0 \o)} = 0. (7) 

Leading approximation (5)-(6) generates the linear iteration scheme 

G = G< > + + • • • + G (n) + • • • , 

where n-th step functional G^> is a solution of the equation 

G W+ t 0— —+»(/{ — tr— 75 — tr^— — } =V*— X • (8) 

or] ^dr] orj or] or] > orj 

A solution of Eq.(8) is the functional 

Q(n) _ p(n)Q(0) 

where is a polynomial of 2n-th degree on source rj. 

It follows from (6), that the unique connected Green function of the leading ap- 
proximation is the propagator S^°\ Other connected Green functions appear in the 
following iteration steps. 

A solution of equation (7) is the free propagator 

m — p 

with dynamical mass m, which is a solution of the gap equation of the NJL model: 

m = -Aigm [ — ^ . (9) 
J m z — p z 

Here and below dp = d 4 p/(2-7r) 4 . 

The divergent integral in the r.h.s. of equation (9) should be dealt with by some 
regularization. The chiral-symmetric, trivial solution m = 0, always exists. Physically 
preferable (energetically favorable) solution is a solution with which corresponds 

to DCSB. Below we shall consider just such a solution. 

The first iteration is the functional 

G« = {^Tr(^ 1} *r? 2 ) +Tr(s« *i/)}g<°>. 

Taking into account Eqs.(5)-(8), we obtain the first iteration equation for the two- 
particle function 

S? ( l I ) = S^(x - y>)SM(x' - y)+ (10) 

+ l gJdx l {{S^{x-x l )S^{x 1 -y))trS < i ) ^ l x j 
-(S^(x - x 1 )75^°)(x 1 - y))tv l5 S? ( X x ) X y ))} 
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and the equation for the first iteration corrections to the propagator 

s », ( , -,>=* / „ lS «»>(, - „)K> ( i\ » ) - ( » ) , 5 } + (id 

+i 5 /" dx 1 S {0) (x-x 1 )S io) (x l -y)trS W (0). 
A solution of the linear integral equation (10) is 

s?( x , y ,\ = 

\x' y' ) 

= [ dx 1 dy 1 dx' 1 dy' 1 S^(x-x 1 )S^(x'-x' 1 )F 2 ( X ) y ) ) S ^ ( Vl -y)S^ (y'.-y'), (12) 
J \ x i Vi J 

where the amputated two-particle function F 2 is 

f 2 (* y y ^=-[s ( % 1 ( X -y f )[s^]- 1 ^-y)+ (13) 

+5(x - y)5{x' - y'){\ ® 1 • A a {x - x') + 75 <g> 75 • A n (x - a/)}, 
and scalar amplitudes A a and A n in momentum space are 

Mp) = ~ 1 _ 19 l y L s (p)=ig JdqtiS^(p + q)S ( -°\q), (14) 



and 



AAP) = l + l [ p{p y L P(P)=*9 J dq tr (p + q) l5 S^ (q) l5 . (15) 

Taking into account the gap equation (9) with m/0, one can obtain (in translational- 
invariant regularization) for A a and A^ the following representations: 

2 (4m 2 — p )Io(P) 

AM = WWY (17> 

Here ^ 

hip) = / — - 2 -. (18) 

7 (m z — (p + q) z ){m z — q z ) 

Eq.(ll) for taking into account the results for S^, reduces in the momentum 
space to a system of simple algebraic equations. Introducing the first iteration mass 
operator = [S^]- 1 ★S'W * [S"^]- 1 we obtain from equation (11): 

^(x) = ig5(x)tTS^(0) + S^(x)A f7 (x) + S (0 \-x)A 7r (x). (19) 

Let us briefly discuss the following step of the iteration scheme. A solution of the 
second iteration functional-derivative equation has the form 

G {2) = {^(sf * , 4 ) + i T r(sf * ^) + i T r(sf * , 2 ) + Tr(s<( 2 ) * r,) } G (°>, 

(2) 

i.e., in the second iteration equations a four-particle and three-particle functions (S4 
and S^) appear. Equations for and S^ 2 -* have the same form as for the first 
iteration, except for inhomogeneous terms, which contain S\ and S3 for the second 
iteration equations. 
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3 Dimensional regularization in the NJL model 



Due to non-renormalizability of the NJL model, the regularization is an essential com- 
ponent of the model. 

We shall use in this work the dimensional regularization in a variant form, which 
was proposed in [11]. The general principles of the approach of [11] are the following: 

• All calculations are made in 4-dimensional Euclidean space; 

• Translational invariance is assumed; 

• The regularization procedure consists in transformation of integration measure 
by a weight function which provides a convergence of the integrals. 

In this approach the dimensional regularization is, in essence, a variant of analytical 
regularization. This point is very important for the use and interpretation of this 
regularization. In this connection we shall, in what follows, use the term "analytical" 
for this regularization to stress its peculiarities in contrast to the usual treatment of 
the dimensional regularization in a formal D-dimensional space. Let us consider the 
gap equation (9) of the NJL model as an example. 

The gap equation with m / in Euclidean space after angle integration is 

^4 f qldq 2 e 



y 2tt YJ 




(2vr) 4 y m 2 +q 2 ' 

where = 2ir 2 is a surface of 4-dimensional unit sphere. In accordance with the 
previous comment, we modify the integrand by the weight function: 

WAMle) = w A (q 2 e )w D (q 2 e ) = 6{A 2 - 

The weight function w\ t o is the product of two weight functions w\ and wd- The func- 
tion wa corresponds to the four dimensional cutoff regularization, while the function 
wd corresponds to dimensional (analytical) regularization. 
A calculation of the integral over dq 2 results in 

Here B x (u,v) is the incomplete Beta- function, 
a) Cutoff: By setting D = 4, we have 



m 2 A 2 \ 



where ka = 9 A 2 /Air 2 . 

This relation coincides exactly with the classical result of [1]. 

b) Analytical (dimensional) regularization: When A 2 — ► 00, taking into account the 
formula 

Bi(D/2, 1 - D/2) = r(D/2)r(l - D/2) 
and rescaling the scale parameter /i 2 as 

((i2)2 . D/2 _ f^Mf-or, m 
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we obtain for the gap equation 

/ o \ D/2-2 

(21) 




4ttM 2 J 



Here instead of g we introduce the dimensionless quantity 

Equation (21) corresponds exactly to the calculation of the initial integral with the 
formal prescription of D-dimensional integration 

,~ = J*q_ (M 2 ) 2 ' D / 2 d D q 
q ~ (2vr) 4 ^ (2tt) d 

but in our case the calculation was performed in the usual 4-dimensional space, i.e. in 
our treatment D is not a dimension of some space, but a parameter which provides the 
convergence. In particular, we are not constrained by the limit D — > 4 for the analysis 
of the results. We suppose that a possible treatment of regularization parameter is a 
power of some additional factor, which is a measure of gluon influence on the effective 
local four-quark self-interaction of the NJL model. 

Below we shall use the regularization parameter £ as 2 

D = 2 - 2£. (23) 

In terms of the parameter £ gap equation (21) has the form 

The region of convergence of the integral is < £ < 1. As we shall see (Section 5), 
this is also the region for the physical values of the model parameters. 
Chiral condensate in the leading approximation is 



(c 3 )(°> = itr 



4vr 2 ^ 



m? 




Integral 1$ (see (18)), which is a part of scalar amplitudes A a and A n , can also 
be calculated as above. Transforming to Euclidean metric, introducing a standard 
Feynman parameterization, and changing an integration variable (which is possible 
due to translational invariance of the procedure, see [11]) we can perform the angular 
integration. According to the procedure, we introduce into the integrand, a weight 
function wd{cl 2 ) and, after the same rescaling (20), obtain the result, which also exactly 
corresponds to the result of integration with the formal transition to D-dimensional 
space. 

Taking into account the gap equation (24), we obtain the pole approximation for 
the scalar amplitudes in analytical regularization 

A (v) ~ 1 = 2fr™ 2 (i + 2fl 

a{P) - 2(4m 2 -p 2 )/ (4m 2 ) (4m 2 - p 2 )^ ' 1 ] 



2 Note that parameter £ is different from the commonly used parameter e = (4 — D)/2. They are connected 
by the relation e = 1 + £. Introduction of this notation prevents unnecessary associations with the usual 
treatment of dimensional regularization. Furthermore, in terms of the parameter £ all subsequent formulae 
of the NJL model acquire simple forms. 
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1 



2igm 2 



(26) 



2p2/ ( ) > 

which correspond to contributions of sigma-meson and pion. We shall use these ex- 
pressions in the following calculations. 



4 Meson contributions to the chiral condensate 



As a measure of quantum fluctuations of the chiral field, consider a ratio of first iteration 
condensate 

(c 3 ) {1) =itTS w (0) (27) 



to the leading-approximation condensate (c 3 )^: 



( c 3)(l) 

r = 



( c 3)(0) 



(28) 



Here r a is a contribution of the scalar meson (sigma-meson) and r w is a contribution 
of the pseudoscalar meson (pion). 
From Eq.(19), we obtain 



r„ = — 



2ig r dpdq(m 2 + 3p 2 - 2{pq)) 
£ J {m 2 — p 2 ) 2 {m 2 — {p — q) 2 ) 



and 



dpdq(m 2 — p 2 + 2(pg)) 



-^r(g). 



£ J (m 2 — p 2 ) 2 (m 2 — (p — q) 2 ) 

Calculation of these two-loop integrals gives us the desired answer. 

The integrals for tv are calculated in the analytical regularization in closed form 
and give us a very simple expression 



r ' = TC 

The scalar contribution can be represented as 



J o 



du 



4*r(§ + p 

20rT(3 + Oio (4-3u) 1 +€ 



(3-2n)(l-Oi ? (l + e,2-^;3 + e 



(n - 2) 2 
4-3u 



(29) 



)" (30) 



-(l + 20F(l+C,l-e;3 + ^^-^) 



where F(a, b; c; z) is the Gauss hypergeometric function. 

Note, that r n and r a depend only on the parameter £ and do not depend on other 
parameters of the model. This feature is a peculiarity of the analytical regularization. 

Results of the calculations in the region < £ < 1 are shown in Table 1. One can 
see that the contribution of the sigma-meson is small, while the contribution of the 
pion in the region < £ < 0.3 4- 0.5 is significant and, strictly speaking, in this last 
region we cannot consider this contribution as a mere correction, i.e., at such values 
of the parameter ^ the quantum fluctuations are large and can lead to a principal 
modification of the entire physical picture of the NJL model (in spirit of work [15], for 
example) . 
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c 


fa 


r-K 


„(U /JO) 


0.1 


0.264 


2.50 


0.556 


0.2 


0.189 


1.250 


0.346 


0.3 


0.119 


0.833 


0.250 


0.4 


0.057 


0.625 


0.189 


0.5 





0.50 


0.145 


0.6 


-0.050 


0.417 


0.110 


0.7 


-0.094 


0.357 


0.081 


0.8 


-0.131 


0.313 


0.057 


0.9 


-0.162 


0.278 


0.037 


1.0 


-0.188 


0.250 


0.021 



Table 1: Relative contributions of sigma-meson (r a ) and pion (r n ) to the chiral condensate 
c 3 of U(l)-model and a ratio (c^/c^) of first iteration condensate to leading-approximation 
condensate as functions of the regularization parameter £. 



5 SU(2)-model 



The U(l)-model has no direct physical applications, and we cannot estimate the phys- 
ical values of model parameters in the framework of this model and, thus cannot reach 
any definite conclusion about the meson contributions. 

To study the possible values of the parameters m, k and £, let us consider a model 
connected with light hadron phenomenology. Accordingly we study the model defined 
by the Lagrangian 



ipidi/j + H 



(<K0 2 + (^75T a Y>) 5 



(31) 



1,2 



Here ip = ^ < j'°\ a = 1, 2, 3, 4 is Dirac spinor index; c = 1, . . . , n c is color index; j 
is isotopic (flavor) index; r a are generators of SU(2) (Pauli matrices); a = 1, 2, 3. This 
model has the chiral symmetry of SUv(2) x 5C/yi(2)-group. We call this the Sub- 
model. 

Mean-field expansion in bilocal-source formalism is constructed with the same 
scheme as in the U(l)-model above, and so without going to the details, we enumerate 
only the major differences from corresponding results of U(l)-model. 

In the leading approximation, the propagator is diagonal over color and flavor: 



? (0) 
cd,jk 



$cd5jk(m-p) 1 
and the gap equation for the SU(2)-model is 

1 = -Sign c j 



dp 



P 



(32) 



(33) 



The chiral condensate is defined for each flavor, i.e., 

c 3 =< 0\uu\0 >, c 3 d =< 0\dd\0 > . 

(In the chiral limit Cd = c u .) Two-particle first iteration amplitude A (connected part 
of amputated two-particle function S^) has the following color and flavor structure: 



A cd,jk _ srcdrc'd' 
A dd'j'k' — 



(34) 



8 



Scalar amplitudes A a and A n are denned by formulae (14) and (15), but traces in the 
definitions of loops L$ and Lp are now taken over all discrete indices. 

First iteration mass operator T,^ = [S^]^ 1 * * [5^°^] _1 is diagonal over color 
and flavor, and is connected with the scalar amplitudes by relation 



Z (1 \x)f k = S cd 5 jk ■ ig5{x)tvS (1 \o) + S (0 \x)A a (x) + 35 (0) (-x)^(x) 



(35) 



For the ratio of the first iteration condensate to the leading-approximation conden- 
sate we obtain 

r = r a + TV = (36) 

Mgn c f dqdp 
1 — 8i gn c J J (to 2 — p 2 ) 2 {m 2 — (p — q) 2 ) 

where 



/ ( m 2-p2 ) 2ff_ (p _ g) 2) ((3P 2 -2( M )+TO 2 )^( 9 )+3(TO 2 V+2(M))4r(5)) , 



f TO +p z 

J = dp- 



(to 2 — p 2 ) 2 

The gap equation in analytical regularization for the SU(2)-model has exactly same 
form as in (24), if we redefine the dimensionless constant k as 

gn c m 2 

k = ^ 2 = 2n c K (37) 

(Here and henceforward, the index denotes corresponding quantities of the U(l)- 
model.) 

The scalar amplitudes of SU(2)-model can be written as 

a An c (Am 2 — p 2 )Io ' ^ ^ 

A * = 1 ~2~F' ( 39 ) 

An c p A h 

where Iq is defined by formula (18). 

Taking into account the above relations and definitions we obtain for the meson 
contributions of the SU(2)-model: 

13 

r a = - — r 0a , r w = - — r 07r , (40) 
2n c 2n c 

where r§ a and are corresponding contributions of the U(l)-model. 

For physical values of colors n c = 3, both contributions decrease compared to 
the corresponding contributions of the U(l)-model: r a to six times, and r n to half. 
Accordingly the bounds of the region of large fluctuations are moved. Results of the 
calculations in region < £ < 1 are shown in Table 2. We see, that for the Sub- 
model at £ > 0.2 the first iteration condensate is not more than 20% different 
from the leading-approximation condensate, i.e., we may consider these values of the 
regularization parameter to be in the stability zone with respect to fluctuations caused 
by the meson contributions. 
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c 


fa 


r-K 


„(U /JO) 


0.1 


0.044 


1.250 


0.319 


0.2 


0.032 


0.625 


0.183 


0.3 


0.020 


0.417 


0.128 


0.4 


0.010 


0.313 


0.098 


0.5 





0.250 


0.077 


0.6 


-0.008 


0.209 


0.063 


0.7 


-0.016 


0.179 


0.052 


0.8 


-0.022 


0.157 


0.043 


0.9 


-0.027 


0.139 


0.036 


1.0 


-0.031 


0.125 


0.030 



Table 2: Relative contributions of sigma-meson (r a ) and pion (r n ) to chiral condensate 
c 3 of SU(2)-model and ratio of first-step condensate to leading-approximation condensate 
(cW/cW) as functions of regularization parameter £. 



6 Choice of parameters 

In determining the SU(2)-model parameters (dynamical quark mass m, regularization 
parameter £ and coupling g (or dimensionless constant k = gn c m 2 /2tt 2 )), it is necessary 
to connect the values of the parameters with measurable quantities. We accordingly 
choose the values of the pion decay constant, the chiral condensate and the decay width 
of 7r°-meson into two photons: 7r° — > 77. 

Pion decay constant /„•= 93 MeV is defined by 

=< 0^ 7/1 7 5 y ij>\b, k >, (41) 

where \b, k > is a pion state b with momentum 

For the SU(2)-model the following formula for f w (see, for example, [4]): 

fl = -4m c m 2 / (0) (42) 

is regularization- independent. In the analytical regularization we have io(0) = i$,/l6ir 2 K 
and obtain: 

fl = |- (43) 

The leading-approximation chiral condensate is 

c = (< 0|^V|0 > /2) 1/3 = -{m/2g) l / z , (44) 

which is regularization-independent. Decay width r 7r o 77 = 7.7 KeV has been calculated 
in the analytical regularization for SU(2)-modcl in ref.[ll]. In our notation, this formula 
has the form 



64vr 3 /2 K 2 

Here a = 1/137 is the fine structure, m^= 135 MeV is the mass of 7r°-meson. 

Formulae (43)— (45) enable us to determine the SU(2)-model parameters. 

Chiral condensate c = (< tpifj > /2) 1 / 3 is not a directly measurable quantity. We 
determine the parameter values for some typical values of the chiral condensate. 
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At c = -160 MeV we have £ = 1, m = 475 MeV, k = 2. Note, that for this 
value of regularization parameter £, the correction to condensate c is about 3% (see 
Table 2), i.e. the model is stable with respect to quantum fluctuations caused by 
the mesons. On the other hand, such a low value of the chiral condensate hardly 
agrees with phenomenology, since the Gell-Mann-Oakes-Renner formula leads to large 
values of current quark masses. At c = —200 MeV we obtain £ = 0.44, m = 400 
MeV, k = 0.62. The correction to condensate c is 9%. At c = —250 MeV, we obtain 
£ = 0.2, m = 370 MeV, k = 0.24, and the correction to condensate c is more then 
18%. Thus fixing the model parameters with formulae (43)— (45) at phenomenologically 
acceptable condensate values c = —(200 -j- 250) MeV corresponds to the condensate 
corrections of the order of 10 -j- 20%. 

The model parameters as chosen above, were implemented in the leading-approximation 
formula (44) for the chiral condensate. The calculated first iteration correction mod- 
ifies the choice of parameters by the following modification of the chiral-condensate 
formula: 



(Accordingly, (43) and (45) are the same with the substitution m — > m*, g — > g*, k — > 



At c = -200 MeV: £* 0.56, m* 420 MeV, k* 9* 0.86; the condensate correction 
is 7%; 

At c = -250 MeV: £* = 0.3, m* ^ 380 MeV, k* = 0.39; the condensate correction 
is 13%. 



Upon comparison of the given values with the ones above, the modification of the 
parameter choice with formula (46) diminishes the relative variation of the condensate, 
i.e., stabilizes the situation. It is a consequence of the positivity of the principal (pion) 
correction to the leading-approximation condensate. 



According to the results obtained, the analytically regularized NJL model gives us 
simple closed formulae not only for the scalar amplitudes and the pion decay constant, 
but also for the pion contribution to the chiral condensate. As it follows from our 
results, in the analytically regularized NJL model, this pion contribution is significant 
and should be taken into account with a choice of physical values of the model param- 
eters. At the same time, the presented analysis demonstrates that in the analytically 
regularized NJL model, for physical values of the parameters, taking into account next- 
to-leading order of the mean-field expansion, does not lead to any pathologies such as 
a disappearance of the DCSB order parameters or of the Goldstone boson in the spirit 
of ref.[15, 16]. 

Acknowledgment 

Authors are grateful to Profs. S.A. Hadjiev and K.G. Klimenko for useful discussions. 




(46) 
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